Exact relations between multifractal exponents at the Anderson transition 
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Two exact relations between mutlifractal exponents are shown to hold at the critical point of the 
Anderson localization transition. The first relation implies a symmetry of the multifractal spectrum 
linking the multifractal exponents with indices q < 1/2 to those with q > 1/2. The second relation 
connects the wave function multifractality to that of Wigner delay times in a system with a lead 
attached. 
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Anderson localization transitions in disordered elec- 
tronic systems represent a remarkable class of quantum 
phase transitions. Here we understand the term "Ander- 
son transition" in a broad sense, including both the loca- 
lization-delocalization transitions and the quantum Hall 
transitions between two phases with localized states. A 
hallmark of these transitions is the multifractality (MP) 
of electronic wave functions, describing their strong fluc- 
tuations at criticality 0: the wave functions are charac- 
terized by a whole set of fractal dimensions Dq different 
from the spatial dimensionality d. While at present the 
wave function MP is routinely observed only in computer 
simulations, rapidly developing imaging techniques allow 
to hope for its forthcoming experimental observation. 

In this Letter we point out two exact relations satisfied 
by the multifractal dimensions. The first of these rela- 
tions connects exponents with q larger and smaller than 
1/2. The second relation links the multifractal indices 
for the wave functions of a closed system to those for the 
Wigner delay times characterizing the wave scattering 
from the same system via an attached lead. 

We begin by considering a functional relation for 
the distribution function of the local density of states 
(LDOS) p. 



(1) 



Here p is the LDOS normalized to its average value, 
p = p/{p) (the normalization factor is not critical and 
plays no role for our discussion). What is of central im- 
portance here is the status of Eq. Specifically, this 
formula is exact on the level of the non-linear cr-model 
(NLcrM). It was derived for the first time in Ref. 2\ for 
the case of systems with broken time-reversal invariance 
(unitary symmetry class). The reasons for its general 
validity were revealed in Ref. and an explicit deriva- 
tion for all standard Wigner-Dyson symmetry classes was 
provided in Ref. Q . 

Clearly, a mapping of a particular microscopic model of 
a disordered system (e.g., continuous model of an electron 
moving in a random potential with certain correlation 
function, or a discrete tight-binding model, etc.) onto the 



NLcrM is not exact. More specifically, it is approximately 
valid in the case of weak disorder and breaks down for 
strong disorder. Therefore, the relation has the same 
status. Nevertheless, we argue that a relation between 
the anomalous multifractal exponents characterizing 
the behavior of the moments (p') at criticality. 
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which follows from (Q, is exact. 

Before explaining this, we digress with a brief reminder 
of the wavefunction MP formalism; the reader is referred 
to the reviews Q for a more detailed exposition. The 
moments of a wave function (so-called inverse partici- 
pation ratios) Pq ~ J d''r|'0(r)p'^ show at criticality an 
anomalous scaling with respect to the system size L, 



Tq = Dqiq-l)^diq-l) + l 



(3) 
(4) 



Here Ag are anomalous exponents distinguishing a crit- 
ical point from a metallic phase [where Tq = d{q — 1)]. 
These exponents also govern the scaling of the moments 
of LDOS, 



(5) 



Equivalently, the MP can be described by so-called sin- 
gularity spectrum /(a), which is the Legendre transform 
of Tq. Its meaning is as follows: the average measure of a 
set of those points r in a sample, where the wave function 
behaves as |^/'^(r)| ^ L^", scales with L as L-^^^K 

We return now to the proof of exactness of Eq. |(5J. 
The central argument relies crucially on the universality 
of critical properties at the Anderson transition. Specif- 
ically, while the mapping of the original microscopic 
model onto the NLtrM is at most approximate, one can 
find another microscopic model (e.g. A^-orbital Wegner 
model in the limit TV ^ oo , which can also be viewed 
as a model of a granular metal 0) that can be reduced 
exactly to the NLcrM. In view of the universality, the orig- 
inal microscopic model and the NLcrM must flow into the 
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same fixed point in the infrared limit and will thus have 
the same critical exponents. Therefore, the relation ^ 
must hold not only in the NLcrM approximation, but be 
an intrinsic property of any generic microscopic model, 
even though the validity of Eq. is in general only ap- 
proximate. 

The moments of the LDOS and of the wave function 
intensity, which we considered above, are properties of a 
closed system. An alternative method to study the local 
properties is to open the system by attaching a perfectly 
coupled single-channel lead at a point r. The system 
can then be characterized by the Wigner delay time tw 
(energy derivative of the scattering phase shift), whose 
statistical properties attracted a lot of research interest 
in recent years, see 0, |^ ■ For convenience, we will con- 
sider below the dimensionless delay time tw — iwA/27r 
normalized to the mean level spacing A. At the Anderson 
transition point the corresponding distribution function, 
Vw{iw), will reflect the criticality of the system 

To establish a connection between the MF of wave 
functions and that of delay times, we recall a relation 
between Vw and the distribution function Vy of normal- 
ized wave function intensities y = l/|?/'^(r)| {V ^ L'^ is 
the system volume), 

Vw{iw)=i^Vy{i^)- (6) 

This formula was derived in Ref. and has the same 
status as Eq. it is exact on the level of the NLcrM. In 
particular, it irnplies the corresponding relation between 
the exponents @| 

Iq = n+g, (7) 

where the indices 7g characterize the scaling of moments 
of the inverse delay time, {i^) oc L^^i . Applying 
the same argumentation as used above for derivation of 
Eq. we conclude that the relation Q must again be 
exact for any generic microscopic model. 

The following point should be emphasized here. 
Strictly speaking, the moments (i^*) with q < —3/2 are 
divergent for the one-channel scattering problem. To de- 
fine the exponent 7g for this case one should consider a 
lead with several conducting channels. This is analogous 
to the coarse-graining procedure for defining the wave- 
function exponent Tg with negative q discussed below. 
Equation (0 holds for such negative q < —3/2 as well, 
by analytical continuation. 

We turn now to the analysis of consequences and ap- 
plications of the derived relations, mainly concentrating 
on Eq. |5J. First, we rewrite the relation ^ in terms of 
the exponents r^, 

Tg - Tl_g = d{2q - 1). (8) 

Second, performing the Legendre transformation, 
f{aq) — qaq — Tg with aq = dTq/dq, we get 

/(2d-a) =/(a)+d-a. (9) 



Equation ((SJ maps the part of the singularity spectrum 
with a < d to that with a > d. A particular conse- 
quence of this is that the support of the singularity spec- 
trum /(a) (i.e. the region where it is different from — oo) 
is bounded by the interval [0, 2d]. The lower boundary, 
a > 0, is a trivial consequence of the wave function nor- 
malization; the upper boundary, a < 2d, follows then 
from our relation Q. 

It is worth mentioning that the results for the /(a) 
spectrum, as obtained numerically for the 3d Anderson 
transition in a number of publications [lol llll| . are in 
contradiction with this upper boundary [and thus with 
the relation 0]. We believe that this is a consequence of 
an incomplete analysis of numerical data in [i(tI . In- 
deed, it was shown recently ^3 that the earlier numerics 
on the wave function MF suffered strongly from the ab- 
sence of ensemble averaging and from finite-size effects. 
The problems become even more severe for negative mo- 
ments, q < 0, corresponding to the large-a part of the 
singularity spectrum. This is evident, in particular, from 
Fig. 6 of Ref. |0| where a strong drift of large-a part 
of /(a) (towards our upper boundary a < 2d = 6) with 
increasing system size is seen. 

Let us analyze the implication of our relation for the 
weak-coupling expansion of the critical exponents that 
can be developed in 2-1- e dimensions (where MF is weak). 
Since Eq. (|2Jl is exact, it should hold in all orders of the 
e-expansion. The known results for the e-expansion of 
Ag up to 4-loop order [3| do satisfy this property. In 
particular, the result for the orthogonal symmetry class 
reads 

Ag = q{l^q)e+^q(q-lMq-l) + l]e^ + 0{e^). (10) 

It is indeed seen that Ag depends on q via the combina- 
tion q{l — q) only, in agreement with the relation 101 . 

As a further application of Eq. we consider the 
model of power-law random banded matrices (PRBM), 
{\H^j\'^) = (l + |i-j|Vfe^)"^ This model (that describes 
a Id system with long-range 1/r random hopping) defines 
a family of critical theories parametrized by < 6 < (X) 
and allows to study the evolution of the critical system 
from the weak- to the strong-MF regime with decreasing 
& llllli- While for & > 1 (weak MF) the PRBM model 
can be approximately mapped to the NLcrM, for small 
b (strong MF) this mapping is not applicable, and the 
multifractal spectrum was analyzed in |l5| by a different 
method. Our statement about the exactness of Eqs. 
remains valid for the PRBM model. Indeed, we can 
construct a "granular" generalization of the model with 
N ^ 1 states at each site of the Id lattice and with 
hopping matrix elements between all states decaying with 
distance r as {b/N)r~^. Changing the overall prefactor 
b in the hopping amplitude will yield a family of critical 
models that should flow in the infrared limit to the same 
line of critical points as the family of PRBM models. In 
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this way, the PRBM model with an arbitrary value of b 
can be associated with an N-orbital model with some b 
that will have the same critical properties. On the other 
hand, the latter model can be mapped onto the NLcrM, 
which allows us again to derive the relations l(2Jl, Q) for 
the critical exponents. 

We have verified the validity of the relation (jj)) by a nu- 
merical simulation of the PRBM model. The exponents 
Tq were extracted from the scaling of the inverse partic- 
ipation ratios (Pq) for system sizes L in the range from 
512 to 4096. The number of disorder realizations was 
ranging from 2 x 10^ for L = 512 to 1000 for L = 4096. 
It should be stressed that evaluation of negative moments 
requires special care, since the inverse participation ra- 
tio, as defined in Eq. Q, is divergent because of zeros of 
the wave function. These zeros, related to oscillations of 
the wave function on the scale of the wave length, have 
nothing to do with multifractal properties characterizing 
smooth envelopes of wave functions. To find Tq with neg- 
ative g, we have first smoothed by averaging over 
blocks of the size m — 16, and then applied Eq. Q. This 
makes finite-size effects (and thus a numerical inaccuracy 
in evaluation of Tq) for negative moments considerably 
more pronounced than for g > 0. 

The results of the numerical simulations for the PRBM 
ensemble with several values of 6, spanning the whole in- 
terval from the weak-MF to strong-MF regime, are shown 
in Figs. ^ 13 The data in Fig. ^nicely confirm the sym- 
metry relation (01. A small difference between A, and 
Ai_g can be considered as a measure of the numerical 
accuracy of evaluation of the exponents. As discussed 
above, the errors are mainly due to moments with nega- 
tive q. In Fig. 121 the same numerical data are presented 
in the form of the singularity spectrum f{a). To demon- 
strate that the data support very well the relation 0, 
we also show the function /(2 — a) + a — 1. 

We will now demonstrate the high utility of Eq.ljSJ 
by applying it for the analytical evaluation of exponents 
with q < 1/2 in the "non-NLcrM" limit, 6 < 1. As was 
found in (1!t| . the multifractal exponents in this regime 
are given for g > 1/2 by 
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2bT{q) , 

2 r(g-i/2) 
r{q-i) 

I ^7T{q-l/2) ' 



(11) 



q^l/2, 
g > 1. 



(12) 



In terms of the singularity spectrum /(a), this means 

/(a) ~ 2bFia/2b), (13) 

where F{A) is the Legendre transform of T{q) with the 
asymptotics 
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-l/irA , 
A/2 , 



A^O , 
A^ oo 



(14) 
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FIG. 1: Multifractal exponents Ag for the PRBM model with 
6 = 4, 1, 0.3, 0.1. The symmetry Q with respect to the point 
q = 1/2 is evident. A small diflference between A, (full line) 
and Ai_5 (dashed) is due to numerical errors. 




FIG. 2: The data of Fig.^in terms of the singularity spectrum 
/(q). Dashed lines represent /(2 — a) + a — 1, demonstrating 
the validity of Eq. Q . 

Equation ((TTJ, lEl) was derived in |l5l | by a real-space 
renormalization-group method valid for q > 1/2. The re- 
lation |j2Jl allows us now to find the multifractality spec- 
trum for q < 1/2. When translated to /(Q!)-language, 
this yields the singularity spectrum for a > 1, 



/(«) 



fi2-a) + a-l 



26 



a/2 , 



1 



7r(2-Q) ' 



a > 26 
26 



(15) 



In Fig. Owe show the MF spectrum of the PRBM model 
for 6 ~ 0.1. The dashed curve yields the a < 1 behavior, 
Eq.(^, while the full fine is the a > 1 result, Eq. l(T5|l . 

In the limit 6 — > the MF reaches its extreme form 
(for the PRBM model the effective spatial dimensionality 



4 




FIG. 3: Singularity spectrum for the PRBM model with 6 = 
0.1. Dashed line: a < 1 behavior, Ea. !]!.'-!^ : full line: a > 
1 result, Eq. 115L following from the relation ll^J; circles: 
numerical data. Some mismatch between the slopes of the 
two curves at a = 1 is related to the fact that the formula 
is valid to the leading order in 6 <^ 1. 

d = 1; we keep d below for generality) 



, 

d{2q-l) 



9 > 1/2 , 
9 < 1/2 , 



or, in terms of the singularity spectrum 
a/2 , 



-oo , 



0<a<2d 
otherwise . 



(16) 
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The following remark is in order here. The earlier anal- 
ysis 13, 01 of the statistics of critical wave functions on 
the Bethe lattice and in large dimensionality d allows us 
to conjecture that in the limit d — > oo the multifractal 
spectrum at the Anderson transition acquires the same 
extreme form (fTBjl . ((T7|l . We stress, however, that this is 
only a hypothesis waiting for a more rigorous verification. 

The second relation we claim to be exact, Eq.Q, is also 
supported by numerical results obtained for the PRBM 
model. Specifically, the numerical data for the scaling 
of the delay time moments confirm (in combination with 
the results of Ref. 0| on the wave function MF) the 
validity of Eq. Q even in the small-6 limit where the 
mapping of the PRBM to the NLctM fails. 

As a final remark, we note that the notion of MF was 
recently extended to the surface of a critical system ^3 • 
While boundary multifractal exponents are different from 
their bulk counterparts, the relations (jSjl and (jZj) remain 
valid also for surface MF. Indeed, it is not difficult to 
check that the derivation of the relations for the distri- 
bution functions (that served as starting points for our 
analysis), (Q) and retain its validity independently 
on the position (in the bulk or near the boundary) of the 
observation point r. The MF of delay times for a lead 
attached to the boundary has in fact been studied nu- 
merically in the PRBM model in Ref. 0; an analysis of 



the surface MF of wave functions and the verification of 
the relation j?)) at the boundary of this system will be 
presented elsewhere. 

To summarize, we have demonstrated that two exact 
relations, Eq.® [that can be equivalently represented in 
the form © or ^] and Eq.®, hold for multifractal ex- 
ponents at the critical point of the Anderson transition. 
We have applied the first of these relations to the mul- 
tifractality spectrum of the PRBM model and verified 
its validity by numerical simulations. A further analysis 
of implications of these relations is of considerable in- 
terest. Another direction of future research is to study 
whether these relations, derived here for three Wigner- 
Dyson classes, have some analogues for unconventional 
(chiral and superconducting) symmetry classes 1^. 

We thank I. Gruzberg for useful discussions. This work 
was supported by the SPP "Quantcn-Hall-Systcmc" and 
the Center for Functional Nanostructures of the DFG 
(ADM, EE), and by EPSRC grant EP/C515056/1 (YE). 



[7: 

[8 

[9: 

[lo; 

[11 

[12 

[13 
[14 

[15 
[16 

[17 
[18 



Also at Petersburg Nuclear Physics Institute, 188350 
St. Petersburg, Russia. 

M. Janssen, Phys. Rep. 295,1 (1998); A.D. Mirlin, Phys. 
Rep. 326, 259 (2000). 

A.D. Mirlin and Y.V. Fyodorov, Phys. Rev. Lett. 72, 526 

(1994); J. de Physique I (France) 4, 655 (1994). 

Y.V. Fyodorov and D.V. Savin, JETP Lett. 80, 725 

(2004) . 

D.V. Savin, Y.V. Fyodorov, and H.-J. Sommers, JETP 
Letters 82, 544 (2005); Y.V. Fyodorov, D.V. Savin, and 
H.-J. Sommers, J. Phys. A: Math. Gen. 38, 10731 (2005). 
F. Wegner, Phys. Rep. 67, 15 (1980). 
K.B. Efetov, Adv. Phys. 32, 53 (1983); K.B. Efetov, Su- 
persymmetry in disorder and chaos (Cambridge Univer- 
sity Press, 1997). 

Y.V. Fyodorov and H.-J. Sommers, J.Math.Phys. 38, 
1918 (1997); T. Kottos, J. Phys. A: Math. Gen. 38, 10761 

(2005) . 

A. Ossipov and Y.V. Fyodorov, Phys. Rev. B 71, 125133 
(2005). 

J. A. Mendez-Bermudez and T. Kottos, Phys. Rev. B 72, 
064108 (2005). 

M. Schreiber and H. Grussbach, Phys. Rev. Lett. 67, 607 
(1991); H. Grussbach and M. Schreiber, Phys. Rev. B 51, 
663 (1995); T. Terao, Phys. Rev. B 56, 975 (1997). 
M. Schreiber, F. Milde, and R.A. Romer, Comp. Phys. 
Commun. 121-122, 517 (1999). 

F. Evers, A. Mildenberger, and A.D. Mirlin, Phys. Rev. 
B 64, 241303(R) (2001). 

F. Wegner, Nucl. Phys B280 [FS18], 210 (1987). 
A.D. Mirlin, Y.V. Fyodorov, F.-M. Dittes, J. Quezada, 
and T.H. Seligman, Phys. Rev. E 54, 3221 (1996). 
A.D. Mirlin and F. Evers, Phys. Rev. B 62, 7920 (2000). 
A. Mildenberger, F. Evers, and A.D. Mirlin, Phys. Rev. 
B 66, 033109 (2002). 

A.R. Subramaniam, LA. Gruzberg, A.W.W. Lud- 
wig, F. Evers, A. Mildenberger, and A.D. Mirlin, 
cond-mat/05 12040 

A. Ahland, M.R. Zirnbauer, Phys. Rev. B 55, 1142 
(1997); M.R. Zirnbauer, J. Math. Phys. 38, 2007 (1997). 



